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Abstract: We investigated the altitudinal variation of Carpinus betulus 
L. in the Hyrcanian forest using leaf macro-morphological and mi¬ 
cro-morphological traits. We collected a total of 1600 leaves from two 
locations. In each location, we sampled six populations along an altitud¬ 
inal gradient ranging from 100 m to 1,150 m. We found that trees in the 
higher elevations have smaller leaf lamina than those in the lower eleva¬ 
tions. In contrast, leaf mass per area was high at low altitudes and in¬ 
creased newly at the higher ones. Stomatal dimension was negatively 
correlated with elevation, while stomatal density was positively corre¬ 
lated with elevation. We also found that two transects showed the same 
plasticity trend. Leaf area showed the highest plasticity, while the number 
of veins showed the lowest plasticity. This study shows that altitude, and 
related temperature and rainfall, represents an important driving force in 
Carpinus betulus leaf morphological variation. Moreover, our results 
suggest that leaf area, leaf mass per area and stomatal density could 
influence the species responses to different ecological conditions. 
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Introduction 

In Iran, the Alborz Mountains, between the Caspian Sea and the 
Iranian Plateau, create an environment extremely propitious to 
the growth of a unique flora in the Hyrcanian Forest. The Hyr¬ 
canian Forest, most closely to forest in a Euro-Siberian zone, is 
rich in broadleaf species (about 80 tree and 50 shrub species), 
while only three to four conifer species. Ecological conditions 
such as soil, temperature, rainfall and solar radiation vary con¬ 
siderably and visibly with increasing altitudes (Geeske et al. 
1994; Lomolino 2001; Austrheim 2002; Wang et al. 2003; 
Korner 2007). In Hyrcanian forest, low-altitude plants have to 
withstand the unfavourable climatic conditions of dry summers, 
with high temperatures, harsh radiation and scant precipitation. 
On the other hand, high altitude plants have to face the adverse 
winter conditions with low temperatures accompained by high 
irradiance. These diverse climatic conditions result in different 
plant morphological adaptations. 

High-altitude species tend to be morphologically and physio¬ 
logically distinct from closely related species from lower alti¬ 
tudes. Altitude has a major effect on leaf morphology and physi¬ 
ology within a species (Korner 1999). In general, leaf morpho¬ 
logical traits, such as leaf length, width and area, are negatively 
associated with increasing altitude (Korner et al. 1986). In con¬ 
trast, leaf thickness increases with increasing altitude (Korner et 
al. 1989; Roderick et al. 2000). Therefore, a study of the rela¬ 
tionships between altitude and plant morphology is extremely 
useful in ecological studies and could be considered by plant 
physiologists, ecologists and palaeo-botanists. On the other hand, 
this high variability can lead to some ambiguities in plant tax¬ 
onomy and makes difficult a comprehensive understanding about 
the possibility of acclimation and adaptation of a species to par¬ 
ticular environmental conditions. Yousefzadeh et al. (2010) re¬ 
ported limited information about leaf variability of Parrotia per- 
sica along an elevational gradient in the Hyrcanian forest, but it 
seems this species was a poor choice because P. persica doesn’t 
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occur across many elevations. Hornbeam ( Carpinis betulus ) is 
one of the most valuable tree species in the Hyrcanian forest and 
occurs across many elevations from the coastal plain at sea level 
to an altitude of around 1800 m in Quercus-Carpinetum and 
Parrotia-Carpinetum communities (Sagheb Talebi et al. 2003). 
Such a wide range of ecological conditions could foster the for¬ 
mation of different ecotypes and/or botanical varieties. Sabety 
(2001), for example, reported Carpinus betulus var. betu¬ 
lus Browicz 1972, Carpinus betulus var. carpinizza (Host) Neilr, 
Carpinus betulus var. parva Radde-Fomin, and Carpinus betu¬ 
lus var. typic Medo, all based upon leaf morphology in the Hyr¬ 
canian forest. 

The primary objective of our study was to evaluate leaf macro- 
and micro-morphological variability of C. betulus along eleva- 
tional gradients and to determine leaf variability in the species. 


Material and methods 

Study site 

The current research was conducted in the Hyrcanian forest, 
northern Iran (Fig. 1). The forest is a wet region; mean annual 
precipitation varies from 2000 mm in western sites to 600 mm in 
eastern ones. Mean annual temperature varies from 15°C in 
western sites to 17.5°C in eastern ones. More than 90% of the 
region is covered by forest brown soils. Generally, rainfall in the 
mountain forests is positively associated with, and temperature 
negatively associated with, increasing elevation (Sagheb Talebi 
et al. 2003). 



Fig 1. Location of the Hyrcanian forest on a map of Iran (left) and 
location of the elevational gradient in the Hyrcanian forest (right). 

Sampling 

Natural populations of C. betulus located in two different loca¬ 
tions of the Hyrcanian forest were sampled for this study. In each 
location, we sampled six populations along an altitudinal gradi¬ 
ent ranging from 100 m to 1,150 m. We collected 20 fully ex¬ 
panded leaves from the outer (light exposed) part of the crown of 
80 trees, resulting in a total collection of 1,600 leaves. Then the 
leaves of each tree were mixed and 10 leaves were randomly 
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sampled for analysis. Only leaves lacking signs of abnormal 
growth, mechanical damage, or pathogen or insect infestation 
were used. 

Macro morphological traits 

Eleven leaf macro-morphological traits were used to evaluate the 
level of leaf variability among the populations. All of the studied 
leaf traits and their abbreviations are explained in Fig 2. We also 
evaluated four leaf-character ratios: leaf length/leaf width 
(LL/LW), leaf length/petiole length (LL/PL), leaf length/distance 
from leaf base to the widest point (LL/BW), distance from leaf 
base to the leaf widest point/petiole length (BW/PL). These ra¬ 
tios, which form independent shape variables, have been used 
extensively in leaf morphometries (Dichinson et al. 1987; Zaraf- 
shar et al. 2010). Leaf area (A, cm2) was also measured with a 
planimeter. Leaves were then dried at 65°C and weighted on an 
electric balance for the calculation of leaf mass per area (LMA, 
mg cm' 2 ). 



Fig. 2 Presentation of the assessed leaf morphological variables. 

Lamina length (LL), lamina width measured at the base (LW 0.1) and at 
the top (LW 0.9), petiole length (PL), distance from leaf base to the leaf 
maximum width (BW). Number of leaf teeth (T) in bold numbers and 
number of leaf veins (V) in non bold numbers. DBV is the distance 
between the first and the second veins. 

Micro-morphological traits 

In the leaf micro-morphological study, stomatal traits: length 
(LS), width (WS), area (AS) and density (SD = number/mm 2 ), 
were measured. For assessing stomatal traits, five leaves were 
sampled from each tree. Leaves were boiled for ten minutes, then 
very thin sub-samples from the bottom surface of the leaves were 
prepared, with a cutter (Blade). To remove the chlorophyllus 
effects, the thin samples were soaked in sodium hypochlorite. 
The epidermal samples were assessed by light microscopy (LM) 
and photos of stomata were taken with 40X lens. Image Tools 
software (version 2. 0) was used for stomata morphometries. 
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Statistical analysis 

After a normality test, the data were log transformed (base of 
natural logarithms). To generate uncorrelated morphological axis, 
we performed a Principal Components Analysis (PCA) from the 
matrix of Loge-transformed tree means. The aim of generating 
uncorrelated axis was so that any covariance between ecology 
and morphology along one morphological axis could be inter¬ 
preted separately from co-variation along another axis (Bruschi 
2010). Morphological responses to both regions and elevation 
were analyzed by performing General Linear Model (GLM), 
two-way analysis of variance (ANOVA), on the PCA scores with 
transects by locations as the main effects of the analysis. Then 
differences among means were determined by Tukey’s test. Ad¬ 
ditionally, a stepwise linear discriminant analysis was employed 
on all data (using the method of inclusion and removal at each 
step) to determine the relative importance of characters as dis¬ 
criminators and to assess the degree of separation of different 
groups by multivariate measurements. A Pearson correlation 
analysis was used to identify whether the altitudinal gradient of 
temperature affected PCA scores. All statistical analyses were 


conducted using the Statistica 6 software package (StatSoft 

2001 ). 

Leaf plasticity 

Total within-transect plasticity (PI) was calculated for each 
parameter using the smallest and greatest mean values PI = 1 - 
(x/X), where x is the smallest value and X is the largest value for 
any given leaf measure (Ashton et al. 1998; Bruschi et al. 2003). 

Results 

Statistical parameters of leaf traits are presented in Table 1. The 
mean comparison among the localities clearly shows that trees in 
the higher elevations have smaller leaf lamina (LL, LW, BW, 
LW0.1, LW0.9 and A) than those in the lower elevations. In con¬ 
trast, leaf mass per area (LMA) was high at low altitudes and 
increased newly at the higher ones. Stomatal dimension was 
negatively correlated with elevation, while stomatal density was 
positively correlated with elevation. 


Table 1. Mean± standard error of the leaf traits for each locality (bold number belongs to the first transect and non bold belongs to the second 
transect) 

Locality _ 

1 (100 m) 

2 (350 m) 

3 (550 m) 

4 (750 m) 

5 (950 m) 

6 (1150 m) 


2 (350 m) 

3 (550 m) 

4 (750 m) 

5 (950 m) 

6 (1150 m) 


LL 

LW 

PL 

BW 

LW1 

LW9 

V 

DBV 

T 

LMA 

8.74 ±0 .37, 

4.57± 0.21, 

1.25dz0.08, 

4.23±0.24, 

2.385± 0.1, 

0.77±0.06, 

11.78±0.38, 

0.54±0.12, 

5.69±0.21, 

9.44±6.66, 

10.29±0.76 

5.12±0.31 

1.48±0.13 

4.54±0.42 

2.84dz.27 

1.07±0.65 

14.66±0.69 

0.56±.52 

5.46±0.29 

10.94ztl 1.74 

9.24±0 .35, 

4.89± 0.19, 

1.62±0.13, 

4.51± 0.22, 

2.71±0.13, 

0.76±0.05, 

12.22zt0.36, 

0.59±0.11, 

6.35±0.20, 

9.37±1.37, 

11.59±0.76 

5.75±0.31 

1.34±0.07 

5.32±0.42 

3.4±.27 

0.75±0.06 

14.66±0.69 

0.58±0.52 

6±0.21 

ll.31dzl.88 

7.72±0 .45, 

3.98±0.25, 

1.19±0.09, 

3.85±0.293, 

2.38±0 .17, 

0.64±0 .07, 

12.2zt0.46, 

0.86±0.15, 

6.04±0.26, 

8.24ztl.44, 

9.53±0.54 

4.65±0.25 

1.49±0.97 

4.63±0.3 

2.54dz 0.19 

0.78±0.46 

13.23±0.49 

0.47±0.37 

6.7±0.21 

9.24ztl.26 

9.31±0 .45, 

4.42±0.25, 

1.16±0.09, 

4.44±0.29, 

2.47± 0.17, 

0.71±0.07, 

12.4dz0.46, 

0.51±0.15, 

5.6±0.26, 

7.22ztl.29, 

8.87±0.59 

4.33±0.24 

1.44±0.1 

4.15±0.32 

2.63±0.21 

0.66±0.5 

13.4zt0.53 

0.44±0.04 

7.08±0.23 

8.49±0.79 

9.09±0 .42, 

4.18±0.23, 

1.27±0.09, 

4.14±0.27, 

2.4±0.16, 

0.56± 0.06, 

14.0±0.43, 

0.46±0.14, 

6.0±0.24, 

9.06±2.10, 

10.49±0.59 

4.61±0.24 

1.27±0.61 

4.89±0.32 

2.39±0.21 

0.75±0.5 

12.4zt.53 

0.46±0.4 

6.35±0.24 

8.99±1.35 

8.17±0.45, 

3.97±0.25, 

1.14±0.09, 

3.92±0.29, 

2.40±0.17, 

0.76±0.07, 

12.40±0.46, 

0.46±0.15, 

5.76±0.26, 

11.89±3.79, 

8.25±0.34 

3.77±0.34 

1.27±0.06 

3.67±0.19 

2.08±0.12 

0.67±0.02 

13.41±0.31 

0.39±0.2 

5.68±0.13 

10.58±2.44 


Locality 

W 

A 

LLLW 

LL/PL 

LL/BW 

BW/PL 

SL 

WS 

AS 

SD 

1 (100 m) 

256.6±27.4, 

357.4zt41.14 

27.9±2.2, 

33.12±2.9 

1.9±0.05, 

2±0.12 

7.08±.47, 

6.97±.75 

2.08±.05, 

2.26±.l 1 

3.41±0.21, 

3.08±0.37 

26.18dzl.7, 

26.13±1.1 

22.63±.91, 

21.31±.31 

472.95±44.41, 

244.43±68.85 

239.73±69.21, 

384.34zt42.1 


232.4±25.9, 
313.4±41.14 
171.5±33.5, 
211.5±29.09 
235.5±33.5, 
249±31.8 
222.2±31.0, 
290.8±31.8 
263.43±33.56, 
184.28±18.39 


24.7±2.1, 

34.21±2.93 

20 . 8 ± 2 . 8 , 

25.28±2.07 

31.6±2.8, 

27.65±2.27 

24.6±2.5, 


1.9±0.05, 

2.05±0.8 

1.9±0.06, 

2.05±0.86 

2.1±0.67, 

2.05±0.94 

2.18±0.06, 


25.04±2.27 2.27±0.94 
21.6±2.8, 2.06±0.06, 
19.14±1.32 2.19±0.05 


7.2±0.4, 

7.25±0.75 

6.6±0.5, 

6.44±0.53 

8.2±0.5, 

7.21±0.58 

7.1±0.5, 

7.68±0.58 

7.1±0.5, 

6.53±0.33 


2.03±.05, 
2.18±.l 1 
2.05±.07, 
2.05±.08 
2.10±.07, 
2.14±.09 
2.2U.06, 
2.15±.09 
2.11±.07, 
2.28±.05 


3.57±0.20, 

3.31±0.37 

3.24±0.26, 

3.13±0.26 

3.92±0.26, 

3.35±0.28 

3.23±0.24, 

3.57±0.28 

3.40±0.26, 

2.9±0.16 


25.18dzl.74, 

32.20±1.26 

26.84ztl.74, 

25.85±1.26 

25.86±1.56, 

25.84ztl.33 

37.49±1.65, 

25.19±1.26 

25.86±1.56, 

23.87±0.82 


21.5±1.08, 

31±1.71 

20.61zt.27, 

21.06±1.1 

20.26±.73, 

21.55±1.44 

22zt.5, 

21.9±.20 

18.62zt.72, 

22.04zt.96 


395.81±44.41, 

241.36±76.98 

484.37±44.41, 

228.71±76.98 

416.75i39.42, 

210.37i76.98 

401.05i39.72, 

246.56i68.85 

446.01i39.7, 

192.91i37.34 


331.17169.21, 

395.73147.15 

260.97169.21, 

363.95147.15 
277.19i61.9, 
379.03i47.15 

152.88161.91, 
396.75i42.18 

284.38161.91, 
429.23i22.8 


Note: LL(Lamina length), LW(leaf width), PL(petiole length), BW(distance from leaf base to the widest point), LW.l(lamina width measured at the base), 
LW.9(lamina width measured at the top), V(number of leaf veins), DBV(the distance between the first and the second veins), T(Number of leaf teeth), LMA(leaf 
mass per area), W(Leaf weight), A(Leaf area), SL(stomatal length), WS(stomata width), AS(stomatal area), SD(stomata density). 


We extracted six principal components explaining about 
84.5% of the total variability (Table 2). PC I explained only one 
fourth of total variance (25.5%), followed by the next five PC’s 


(16.7%, 12.4%, 10.7%, 9.6% and 9.5%, respectively). Most 
leaf-size parameters, such as lamina length (LL), lamina width 
(LW), petiole length (PL) and distance from lamina base to the 
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lamina maximum width (BW) were most strongly correlated 
with PC I, as expected. Some leaf shape parameters like LL/PL 
and BW/PL were correlated with PC II. In turn, leaf area (A) and 
leaf mass per area (LMA) were highly correlated with PC III. 
Only distance between the first and the second veins (DVB) was 
significantly correlated with PC IV. The LL/BW ratio was 
negatively correlated with PC V whereas the only variable 
positively correlated with PC VI was the LL/LW ratio. Two-way 
ANOVA on PCI (correlated with leaf size) and on PC III 
(correlated with leaf weight and area) scores revealed that 
differences existed between transect and location (p < 0.05) 
(Table 3). Transect and location did not influence leaf-shape; 
PC2 (correlated with LL/PL and BW/PL) did not differ among 
transects and locations (Table 4). Post hoc contrasts of 
differences in PC scores did not show an easily interpretable 
pattern (Table 5). 


Table 2. Factor loadings carried out on macro-morphological data 



Factor 1 

Factor 2 

Factor 3 

Factor 4 

Factor 5 

Factor 6 

LL 

0.932 






LW 

0.828 






PL 

0.724 






BW 

0.896 






DBV 




0.763 



LMA 



0.934 




A 



0.852 




LL/LW 






0.873 

LL/PL 


0.963 





LL/BW 





-0.945 


BW/PL 


0.917 





Eigenvalue 

3.818 

2.521 

1.836 

1.601 

1.437 

1.427 

Explained variance 

0.255 

0.167 

0.124 

0.107 

0.096 

0.095 


Note: LL(Lamina length), LW(leaf width), PL(petiole length), BW(distance 
from leaf base to the widest point), DBV(the distance between the first and 
the second veins), LMA(leaf mass per area), A(Leaf area). 

Table 3. Two-way analysis of variance of leaf stomata data 

_SS_ df MS F _P_ 

Factorl transect 1.62288 1 1.622882 2.045032 0.158465 

location 2.75743 5 0.551486 0.694941 0.629515 

transect* location 7.85073 5 1.570146 1.978578 0.096595 

Error 42.85293 54 0.793573 


Table 4. Two-way Analysis of Variance carried out on each factor 
extracted 




SS 

df 

MS 

F 

P 


Transect 

6.36023 

1 

6.360233 

9.800678 

** 

Factorl 

Location 

14.66930 

5 

2.933860 

4.520876 

** 

transect*location 

19.34682 

5 

3.869365 

5.962423 

*** 


Error 

44.77813 

69 

0.648958 




Transect 

4.24166 

1 

4.241656 

4.446139 

ns 


Location 

2.98919 

5 

0.597838 

0.626659 

ns 

Factor2 

transect*location 

5.71742 

5 

1.143483 

1.198609 

ns 


Error 

65.82660 

69 

0.954009 




Transect 

8.13366 

1 

8.133660 

10.74659 

** 


Location 

10.49579 

5 

2.099158 

2.77351 

* 

Factor 3 

transect*location 

7.34739 

5 

1.469478 

1.94155 

ns 


Error 

52.22330 

69 

0.756859 




Table 5. Tukey’s Post- hoc summary for comparison of transect and 
locations in relation to leaf macro and micromorphological traits 

transect location Macro Micro 

morphological morphological 




factor 1 

factor2 

factor 3 

factor 1 

1 

1 

AC 

A 

AB 

AB 

1 

2 

AC 

A 

A 

A 

1 

3 

AB 

A 

AB 

AB 

1 

4 

AC 

A 

B 

AB 

1 

5 

AB 

A 

BC 

AB 

1 

6 

ABC 

A 

AB 

AB 

2 

1 

AB 

A 

AC 

AB 

2 

2 

B 

A 

A 

AB 

2 

3 

AB 

A 

A 

AB 

2 

4 

ABC 

A 

A 

AB 

2 

5 

AB 

A 

A 

AB 

2 

6 

C 

A 

A 

B 

When 

we 

analyzed 

micro-morphological 

data, only one 

principal 

component 

was extracted 

and 

this component 

accounted 

for 

more than 63% of variations (Table 6). This 


component was negatively correlated with stomatal area, large 
length stomata and stomata index and was positively correlated 
with stomatal density. Results of ANOVA did not reveal 
significant differences between transect and localities (p > 0.05). 
Pearson correlations performed on all PC scores showed that 
only PC I, concerning micro-morphology, was significantly 
related to the temperature gradient within each transect (R = 
-0.43, p < 0.05). Thus, it suggests a significant pattern by which 
locations from the lower altitudes (higher average temperatures) 
have bigger stomata and lower stomatal density. 

Table 6. Factor loadings carried out on micro-morphological data 


Variable 

factor 1 

AS 

-0.922 

SD 

0.857 

stomata index 

-0.952 

Eigen value 

2.536 

Explained variance 

0.634 


Note: AS(stomatal area), SD(stomata density). 


Table 7 shows the most significant characters in the first three 
canonical discriminant functions, which accounted for 100% of 
the variations. The first function accounted for 77% of the total 
variations and separated locations 5 and 6 from the others (Fig. 
3). The standardized coefficients of DFA 1 (Discriminant 
Function Analysis) were highest for vein numbers and petiole 
length. Function 2 represented another 22% of the total variance 
and tends to separate the locations within each altitudinal 
transect. This function was related to tooth number. Function 3 
accounted for only 1% of the total variations and was a 
petiole-length-susceptibility function. Overall misclassification 
was relatively high (17%) with values of 50%, 100%, 100%, 
50%, 100%, and 100% for locations 1, 2, 3, 4, 5 and 6, 
respectively. 
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Table 7. Standardized coefficient from discriminant analysis 



DFA 1 

DFA 2 

DFA 3 

V 

1.998 

0.306 

-0.173 

PL 

-1.917 

0.354 

-0.686 

T 

0.158 

-0.989 

-0.338 

Eigenvalue 

6.134 

1.602 

0.187 

Explained variance 

0.770 

0.220 

0.010 


Note: V (number of leaf veins), PL (petiole length), T (Number of leaf 
teeth). 
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Fig. 3 Scatterplot of the discriminant canonical analysis. Discriminant 
scores for the first two functions of the means of considered traits for the 
different localities (1-6) along the two transects are shown on the axis. 

Values of leaf plasticity were relatively low for all of the 
macro-morphological traits (total mean value: 0.51 ± 0.12). 
Lowest and highest plasticity mean values were recorded for 
number of veins (V) (PI = 0.34) and leaf area (A) (PI = 0.69) 
respectively. Stomatal traits were not very plastic (total mean 
value: 0.25 ± 0.09), with the highest value recorded for stomata 
density (about 0.39). It seems that the two transects showed the 
same plasticity patterns (Fig. 4). 



Fig. 4 Plot of plasticity ( PI) of macro- and micro-morphological traits 

for the two analyzed transects. The abbreviations along X-axis refer to 
Table 1 


Discussion 

Leaf morphological traits of most woody plants are affected by 
different abiotic factors along environmental gradients, in par¬ 
ticular, several studies documented strong phenotypic variations 
as adaptive response to altitude (Schoettle and Rochelle 2000; 
Hovenden and Vander Schoor 2003; Premoli et al. 2007; Youse- 
fzadeh et al. 2010; Bresson et al. 2011). Morphological and 
physiological adjustments allow individuals to cope with envi¬ 
ronmental factors associated with elevational gradient: these 
include dwarfness, compact habit and small, densely pubescent 
leaves. Our study showed that leaf size in C. betulus decreased 
significantly (Table 1 and Table 2) with increasing altitude in 
both elevational transects, although this trend was more clear and 
linear in the second one (Table 1). This result is consistent with 
the findings of previous studies (Korner et al. 1986; Zhang and 
Marshall 1995). Yousefzadeh et al. (2010) observed a significant 
effect of altitude on width of lamina, base angle, number of 
paired leaf veins in Parrotia persica sampled in the Hyrcanian 
forest. In the same area, Akbarian et al (2011) reported signifi¬ 
cant and negative correlation values of leaf blade, petiole length 
and stomatal length with increasing elevation on Alnus subcor- 
data. These morphological variations among leaves appear to be 
a consequence of adaptive strategies to abiotic factors along 
elevational gradients, such as soil moisture (Beerling et al. 1996), 
air temperature (Panek and Waring 1995) and atmospheric C02 
concentrations (Ehleringer and Cerling 1995; Marshall and 
Monserud 1996). We also found that LMA was higher at lower 
and higher altitudes but decreased at mid altitudes. The higher 
LMA probably results from higher temperatures and summer 
drought at lower altitudes (Sletvold and Agren 2011) and cold 
temperature at higher altitudes (Kao and Chang 2001; Bresson et 
al. 2011). In a large scale study of the Catalan Q. ilex forests 
(Ogaya and Penuelas 2007), higher LMA was observed in drier 
environments, but also in more sunny areas and with lower air 
temperatures. The effect of air temperature was even greater than 
the effect of water availability, especially in the case of cold 
winter temperature that was considered as the most important 
factor affecting LMA values in this species. On the other hand, it 
has been pointed out that leaves with high LMA have a structure 
that reduces photosynthetic rates but also water losses by tran¬ 
spiration (Hultine and Marshall 2000). 

Our results showed an increase in stomatal density with in¬ 
creasing altitude although this increase was not linear along the 
altitudinal gradients (Table 1). Previous studies relating stomatal 
density to altitude obtained mixed results. Some found that 
stomatal density increases linearly with altitude (Korner and 
Cochrane 1985), that it increases and then decreases again 
(Korner et al. 1989), decreases with altitude (Schoettle and Ro¬ 
chelle 2000) or remains unchanged (Woodward 1986; Holland 
and Richardson 2009). An increase in SD might be associated to 
an improved efficiency in carbon dioxide uptake (McElwain 
2004), whereas a decrease in SD might reflect a conserving wa¬ 
ter strategy (Schoettle and Rochelle 2000). Li et al. (2006) 
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showed that stomatal area, stomatal length and index of Q. aqui- 
folioides decreased with increasing altitude above 2800 m al¬ 
though leaf morphological and physiological responses to alti¬ 
tudinal gradients were non-linear with increasing altitude. Of the 
ecological variables they investigated, increasing drought from 
lower to higher altitudes was considered the main stress force. 
The authors concluded that colder soils could reduce water up¬ 
take of root system and induce water stress. Significant Spear¬ 
man rank correlation was observed between temperature and PC 
performed for stomatal traits, suggesting that mi¬ 
cro-morphological analysis could be a useful tool to analyze 
phenotypic adaptation in this species. Moreover, this result could 
be interesting from the perspective of paleontological reconstruc¬ 
tion of altitude (Kouwenberg et al. 2007). In general, different 
environments that represent different microclimates influence the 
variation of leaf quantitative characters (Valladares et al. 2002). 
For example, Bruschi et al. (2003) and Bruschi (2010) reported 
lower values of leaf plasticity for xeric populations of Quercus 
petrea. Although central and eastern parts of the Hyrcanian forest 
do not experience the same climate, we found similar patterns of 
plasticity for each transect. Number of veins, number of teeth 
and stomatal dimensions showed the least plasticity, while leaf 
area and stomatal density showed the highest plasticity. Some 
studies showed that vein density were less sensitive to changes in 
environmental conditions, suggesting that leaf vein traits are 
conservative (Sack et al. 2008; Zhu et al. 2011). Vein traits are 
considered to reflect the gas and water change characteristics 
between the leaves and the atmosphere, which are strongly af¬ 
fected by climate (Zhu et al. 2011). It seems, however, that each 
species showed a special pattern in leaf plasticity. Zarafshar et al. 
(2010), for example, reported the highest plasticity for leaf peti¬ 
ole length of Castanea sativa and Akbarian et al (2011) reported 
the highest plasticity for leaf-apex length of Alnus subcordata. 
On the other hand, leaf traits that show the greatest plastic re¬ 
sponse are most important for leaf functioning in different envi¬ 
ronments (Markesteijn et al. 2007). If this is the case, leaf area, 
leaf mass per area and stomatal density could influence altitud¬ 
inal acclimation of Carpinus betulus to climate change (Aitkin et 
al. 2008). To understand adaptation of this species to altitude, 
more large-scale investigations are needed, also taking ecophysi- 
ological traits into account. Moreover, common garden experi¬ 
ments will be required to understand the environmental and ge¬ 
netic components of variation (Bresson et al. 2011). 
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